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Abstract 

We present a statistical study of the characteristics of very large amplitude whistler waves 
inside the terrestrial magnetosphere using waveform capture data from the Wind spacecraft as 
an addition of the study hy \Kellogg et a/.H2010| . We observed 244(65) whistler waves using 
electric(magnetic) field data from the Wind spacecraft finding ~40%(~62%) of the waves have 
peak-to-peak amplitudes of >50 mV/m(>0.5 nT). We present an example waveform capture of 
the largest magnetic field amplitude (>8 nT peak-to-peak) whistler wave ever reported in the ra- 
diation belts. The estimated Poynting flux magnitude associated with this wave is >300 /xW/m^, 
roughly four orders of magnitude above previous estimates. Such large Poynting flux values are 
consistent with rapid energization of electrons. The majority of the largest amplitude whistlers 
occur during magnetically active periods (AE > 200 nT). The waves were observed to exhibit 
a broad range of propagation angles with respect to the magnetic field, 0° < dtB < 90°, which 
showed no consistent variation with magnetic latitude. These results are inconsistent with the idea 
that the whistlers are all generated at the equator, propagating along the magnetic field, and that 
the observed obliqueness is due to propagation effects. We also identified three types of electron 
distributions observed simultaneously with the whistler waves including beam-like, beam/flattop, 
and anisotropic distributions. The whistlers exhibited different characteristics depending on the 
observed electron distributions. For instance, the whistlers observed with anisotropic distribu- 
tions in the radiation belts had larger Af/f than the rest of the whistlers. The majority of the 
waveforms observed in our study have f/f„ < 0.5 and are observed primarily in the radiation belts 
simultaneously with anisotropic electron distributions. 

1 Introduction 



Whistler waves are one of the most ubiquitous wave modes in plasmas. They have been observed in 



the magnetosphere 


Burtis and Helliwell\ \969\ Cattell et al. 2008 Cully et al. IQQ'i \Parrot et al. 


\2003\\Russell etal.\ 


1969| Santolik et al. |2003|, in the solar wind | Breneman et al. 2010| Neubauer\ 
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and Musmann 1977 1, upstream of interplanetary shocks | Wilson III et al. 2009) , upstream of plan 
etary bow shocks | Bertucci et al. 2005 [ Hoppe et al. 198 1| , and in cometary foreshocks | Tsurutani 



et al. 1987 1 . Whistler waves are a right-hand polarized electromagnetic mode that can propagate 



along the magnetic field or at highly oblique angles as a quasi-electrostatic mode near the resonance 
cone interacting with both ions and electrons. Whistlers can be driven unstable by electron temper- 



ature anisotropics | //as/ii'mo/o aw/ ffimwra 1981 Kennel and Petscheck 1966 1, electron heat flux 



I Gary et aE] |1994[ , and electron beams | Sauer and Sydora\ 2010[ Tokar et a/.[|1984[ \Z,hang et'aL\ 
993 1. For simplicity, we include chorus, plasmaspheric hiss, and oblique whistler-modes in our use 



of the term whistler wave. Because whistlers interact strongly with energetic particles \ Kennel and 



Petscheck 1966 Lyons et al. 1972 1, it has been well accepted that they play an important role in 



global radiation belt dynamics. Thus, whistler waves have been a topic of extreme interest for over 
40 years in magnetospheric physics. 

Time-averaged spectral intensities have been used to estimate whistler wave amplitudes for over 



30 years | Gurnett et al. 1976 1, but averaging can underestimate the instantaneous wave amplitudes. 
Typical whistler time-averaged amplitudes are ^0.5 mV/m for the electric field {Mer edith et 0!] 
2001 1 and ^0.01-0.1 nT for the magnetic field \ Horne et al. 2003 2005) . The underestimated wave 
amplitudes may mask important wave-particle interactions which require larger amplitude waves. 
Recent observations hy \Santolik et al. 1 2003), Cattell et al. \ 2008 1, Cully et al. 1 2008), and Le Contel 



\et a/.|p009| showed that previous time averaged spectral intensity measurements may often under- 
estimate the wave amplitude of terrestrial whistler waves by more than two orders of magnitude. 
Test particle simulations have found that electrons can be accelerated to MeV energies in relatively 



short times by these very large amplitude whistler waves | Bortnik et al. 2008 Cattell et al. 2008 



Omura et al. 2007) . These observations have raised new questions regarding the energization and 



Ufetime of radiation belt particles. 

During an eight year period, the Wind spacecraft went through a number of petal orbits into 
the terrestrial magnetosphere. We report on whistler wave statistics for 13 of those petal orbits, 
building on the study by Kellogg et al. | 2010) , who used an automated search algorithm to find large 



amplitude whistlers. A search by eye yielded significantly more whistler waves than reported in Kel 



logg et al. 1 2010) . This is the first study to examine the statistics of the relationship between wave 
propagation angle and magnetic latitude for these very large amplitude whistlers, finding no corre- 
lation between propagation angle and magnetic latitude. This is the first observation of very intense 
whistler Poynting flux. We also present the first statistical study of the characteristics of the electron 
distributions over the energy range of a few eV to 30 keV associated with large amplitude whistler 
waves inside the terrestrial magnetosphere using waveform capture data simultaneously with three 
different types of electron distributions. The paper is organized as follows: Section|2]introduces and 
outlines the data sets and analysis techniques. Section [3] describes the observations, and Section]?] 
discusses the results and conclusions of our study. 



2 Data Sets and Analysis 

Waveform captures were obtained from the WindAVAVES instrument \Bougeret et al.\\\995\ , using 
the time domain sampler (TDS) receiver, which provides a waveform capture (herein called TDS 
sample) of 2048 points with timespans ranging from ^17 ms to ^1000 ms, depending on sample 
rate. TDS samples are utilized from the fast (TDSF) and slow (TDSS) TDS receivers. The TDS 
receivers have a low frequency cutoff of roughly 3.3 Hz for TDSS and 130 Hz for TDSF. The TDS 
buffer stores and evaluates waveforms based upon their amplitude, keeping only the largest events. 
TDSF samples are composed of two electric field vectors in the XY-GSE plane, defined as E;. TDSS 
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samples are composed of four vectors, either three electric (E^) and one magnetic (B j) or three mag- 
netic and one electric field components. The TDSS samples are rotated into magnetic field-aligned 
coordinates (FACs), determined with magnetic field measurements from the Wind magnetic field 
instrument (MFI) \ Lepping et al. 1995) , where the subscript ||(_L) refer to parallel (perpendicular) to 
the magnetic field. We define the wave amplitude as |E„ | for electric fields and |B„ | for magnetic 
fields. The Wind/WAVES instrument also contains an onboard time-averaged spectral intensity in- 
strument, the thermal noise receiver (TNR), used to analyze low amplitude signals and determine the 
local plasma density using the plasma line when possible. Note that we only use the TNR receiver 
for determination of the local plasma line, not for whistler wave identification. 

The wave vector, k, and propagation angle with respect to the magnetic field, were de- 
termined using Minimum Variance Analysis (MVA) [ Khrabrov and Sonnerup] [199 81 on bandpass 
filtered TDSS samples with three magnetic field components. The frequency ranges for each band- 
pass filter, determined from spectral analysis, were chosen independently for each TDS sample. We 
required intermediate to minimum eigenvalues of the spectral matrix, A„,/Am,„, to satisfy the condi- 
tion > 10.0 if less than 50 field vectors were used in the analysis \Khmbwv and Sonnerup 1998[ . 

Full An steradian particle distributions for electrons and ions were obtained from the Wind/3DP 
EESA and PESA particle detectors for low energies (<30 keV) and the solid state telescopes (SSTs) 



for high energy (-30-500 keV) electrons (SST Foil) and (-70-6000 keV) protons (SST Open) \Lin 



\et al.\ |T995) . Electron distribution functions, from a few eV to 30 keV, were examined for pitch- 
angle anisotropics as a possible free energy source. Higher energy particle distributions using the 
3DP solid state telescope (SST) detectors were used to identify the periods when Wind entered the 
radiation belts. However, since the detectors were built for solar wind use, they often saturate in the 
radiation belts. Thus, they could not be used to look for potential free energy sources or for corre- 
lations between flux enhancements/losses and strong waves. We use the MFI data, combined with 
PESA data, to define the region we call the outer magnetosphere as the region between the terres- 
trial magnetopause and the radiation belts. The magnetopause was defined by the sharp gradient in 
flow speed, density, B;-GSM, and magnetic field turbulence. We defined the radiation belt region by 
the sharp increase(decrease) in omni-directional flux of the high energy electrons (> 100 keV) and 
protons (> 1 MeV) as Wind approached(departed from) perigee in its petal orbits. To allow for the 
possibility of different sources of free energy and/or damping, outer magnetospheric and radiation 
belt whistlers were examined separately. 

We also examined electron temperature anisotropics (T^^/Tny 7^ 1) for both the Eesa Low (j = c) 
and Eesa High (j = h) separately to compare to the threshold criteria for whistler anisotropy insta- 
bility \ Kennel and Petscheck\ ^96()\. The two EESA detectors have energy ranges (<30 keV) well 
below the typical >100 keV radiation belt particles. Thus, using these detectors limits our ability 
to accurately describe the full particle distributions in the radiation belts. However, we will show 
that many of the observed whistlers are in resonance with electrons below 30 keV. Recall that we 
mentioned that the SST detectors often saturate in the radiation belts, thus we are limited to <30 
keV electron distributions in this study. 



3 Observations 

Figure [1] shows an example of a Wind perigee pass from 1998-11-13/13:00:00 UT to 1998-11- 
14/02:00:00 UT. The top panel shows one minute averaged spectral data from the TNR receiver in 
piWI^/Hz- Overplotted is the upper hybrid frequency (f„;,) line. Due to the large amounts of noise 
in the TNR data, the density for the f„,, line was estimated from PESA High density measurements. 
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However, the PESA moments give an underestimate of density at low L-shells, so our density esti- 
mates are likely too low. It can be seen that the whistlers generally occur between L ^3-5 (locations 
marked by vertical lines in all panels), in the radiation belts. Also note that the lowest energy bins 
of the SST detectors saturate between roughly 16:00 and 20:00 UT on 1998-11-13. This is the most 
obvious in the SST Foil panel for the 40 keV line which drops in flux just after 16:00 UT. This 
saturation limited our ability to use the SST detectors when examining the whistler waves. 

Figure [2]plots two whistler waves observed near the equatorial plane, post midnight near L ^ 4. 
These whistlers were obtained from the TDSS instrument (with three electric and one magnetic field 
sampled at 7.5 kHz); the waveforms shown were filtered (frequency range shown in top panels) to 
remove superposed low and high frequency signals. The L-shell, MLT, GSM latitude, lower hybrid 
frequency (f;;, = -//„/„), and f„ for each corresponding whistler are labeled in the figure. Note that 
the 8 nT peak-to-peak amplitude we report here is only one component of the magnetic field. The 
total magnetic amplitude for both waveforms shown in Figure |2] is larger. The polarization of the 
electric fields is elliptical for both examples. The 1998-1 1-13 example is much more oblique (0^ ~ 
73°) and efliptical {X„,JX,„,„ - 10.7) than the 2000-04-10 example (0^ - 39° and /l„„,/A„„, ~ 1.9). 
Omura et al. 1 2007) gave a relationship for the maximum change in kinetic energy given by: 

5.6 X 10"* 8B , 

where L is the L-shell, ^0^ = {£Ieq - (0)1(0 pe^, £Ieq is the equitorial cyclotron frequency, and 
5B/Bo is the ratio of wave amplitude to background magnetic field. For the waves shown in Figure 
|2] (AKE),,,,, ~ 61 MeV(28 MeV) for the wave on 1998-11-13(2000-04-10). The examples show 
that: (1) large amplitude whistler waves in the radiation belts are bursty; (2) electric fields are in 
excess of two orders of magnitude above previous estimates; (3) the first observations of these very 
large amplitude whistler waves with search coil magnetic fields shows amplitudes exceeding two 
orders of magnitude above previous estimates; and (4) the waves are capable of producing electrons 
with energies greater a MeV. 

To illustrate the occurrence rate of these large amplitude whistlers we will discuss two Wind 
perigee passes in detail. During the perigee passes of 1998-11-13 and 2000-04-10, Wind observed 
some of the largest whistler events ever recorded in the magnetosphere. On 1998-11-13 between 
18:15:08.170 UT and 18:21:17.884 UT, or roughly 6 minutes. Wind observed 1 1 TDSF and 3 TDSS 
samples. All 14 TDS samples had |E,„ | > 100 mV/m and all 3 TDSS samples had |B„, | > 4 nT. 
The large number of whistlers observed during such a short duration and spatial region is consistent 
with STEREO observations \Cattell et a/.| [20081 . On 2000-04-10 between 02:45:39.736 UT and 
04:32:44.807 UT, or roughly 1.75 hours, the spacecraft moved through -3.48 hr(-52.2°) of MLT 
and -0.55 L-sheU. During this interval. Wind observed 34 TDS samples, 27 TDSF and 7 TDSS. Of 
those 34 samples, 31 had |E„, | > 80 mV/m and aU 7 TDSS samples had |B,„ | > 1 nT. 

Since we only have four field components for each TDSS sample, we cannot fully describe the 
wave Poynting flux. However, we can calculate part of two components to estimate the magnitude of 
the Poynting flux for these waves. The 1998-11-13(2000-04-10) whistler in Figure|2]has Poynting 
flux mag nitude of >300 AiW/m2(>30 ^iW/m^), which is roughly four(three) orders of magnitude 
larger than the estimates found hy ^Santoh'k et al. | 2010J. To e stimate the possible impact of the large 
Poynting fluxes, we perform the calculation of SantoUk et al. \ 2010| assuming the same estimates of 
>1 MeV electron fluxes, background densities, and field-aligned column area (of a flux tube). We 
find that it would take roughly 5 ms(50 ms) for the whistlers seen in Figure [2] to deposit the neces- 
sary energy density (—10^^ J/m^) to accelerate plasma sheet electrons to 1 MeV, assuming 100% 
efficiency, in the outer radiation belt. If we now assume a 1% efficiency, we find a time scale of 
0.5 seconds(5.5 seconds) necessary for these whistlers to produce the same effect. These estimates 
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are five to six orders of magnitude shorter than the typical estimates. However, as one can see in 
Figure [2j these whistlers are very bursty and their amplitudes are not sustained for longer than lO's 
of milliseconds. We should also note that the larger amplitude waves we observe will likely not 
interact with electrons in the same quasi-linear fashion as described by SantoUk et al. 1 2010 1. 

Figure |3] shows the relationship between peak whistler wave amplitudes and one minute AE- 
Index for the whistlers. The plots show a slight trend of increasing peak amplitudes with increasing 
AE, consistent with increased substorm activity providing more free energy for whistler growth. A 
cursory examination of LANL geosynchronous low energy (30-300 keV) electron data shows large 
injections on the nightside just prior to the two whistlers observed in Figure [2] consistent with the 
AE-index analysis. 

Figure |4] illustrates the different types of electron distribution functions observed at times close 
to and/or concurrently with the whistlers. Each plot represents the parallel and perpendicular cuts 
of the electron distribution function seen as the solid red and dashed blue lines, respectively. The 
first panel. Type A, shows beam-like and/or flattop features in the parallel cut, called a beam/flattop 
electron distribution function. The second panel. Type B, does not show a flattop feature, but does 
show beam-like features. Finally, the third panel. Type C, is an example of an electron distribution 
function that has T^/T|| > 1.0 for ^400 eV to >30 keV, called an anisotropic electron distribution 
function. The anisotropy seen in Type C extended well beyond 1 keV, while the characteristic fea- 
tures of Types A and B occurred primarily below 1 keV. We will refer to whistlers seen near Types 
A and B as beam/flattop-related whistlers and whistlers seen near Type C as anisotropy-related 
whistlers. Note that due to the differences in sampling times (^3 seconds for electron distribution 
functions versus ^17-250 ms for whistlers) there may be short-lived features in the electron distri- 
bution functions that are not observed. Table [T] shows that most of the whistlers were observed in 
the radiation belts and were associated with anisotropic electron distribution functions. 

Figure [5] summarizes the whistler wave location and wave vector direction statistics. The fig- 
ure shows the twelve Wind orbits in the top four panels projected onto the XY-GSM plane with 
color-coded lines and the date for each line is labeled in the corresponding panel (e.g. red line for 
1998-11-13 petal orbit in top left panel). The plots of the petal orbits illustrate the limited coverage 
of the magnetosphere by Wind. Over plotted on each orbit are color-coded symbols, *'s, where the 
colors green, orange, and purple respectively correspond to the following electron distribution func- 
tions: (1) outer magnetospheric beam/flattop, (2) outer magnetospheric anisotropic, and (3) radiation 
belt anisotropic electron distribution functions. For reference, the cumulative number of each type 
of whistler seen in each panel are labeled for corresponding color-coded *'s {e.g. 5 whistlers seen 
with anisotropic electron distribution functions in the outer magnetosphere in top left panel). 

The bottom part of Figure [s] shows two sets of histograms of GSM-Latitude of Wind, Xqsm, 
on the left and S^g, the whistler wave vector angle, on the right. The histogram for Xqsm includes 
239 TDSF and TDSS samples and shows the corresponding range of latitude during the four petal 
orbits. Since could only be determined for TDSS events with three magnetic field components 
and we required that A„„/A„„„ > 10.0 be satisfied, we determined for only 45 whistlers. Roughly 
95%(65%) of the whistlers were observed within 2 R£(10°) of the XY-GSM plane largely due to 
the Wind spacecraft's orbits. All of the outer magnetospheric whistlers for both electron distribu- 
tion function types are observed between 4.5 < L < 12 and all radiation belt whistlers are observed 
between 3 < L < 7.5. Also, 92% of the radiation belt whistlers and all the beam/flattop-related 
whistlers are observed within ±6 hours of midnight, while only 44% of the anisotropic outer mag- 
netospheric whistlers occurred in this range. Most whistlers observed in the outer magnetosphere 
had wave vectors close to parallel (^60% had 0^ < 30°), while the whistlers in the radiation belts 
had a wider range of wave vectors (0° < O^b < 90°). 

Table[T]shows the statistics of all the identified whistler wave samples inside L — 15. Two hun- 
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dred fifty-seven whistlers were observed inside the magnetosphere; of these, 244 whistlers occurred 
between 3 < L < 12. The top part of the table separates the whistlers into three columns by region 
and into two rows defined by the type of electron distribution function. The middle part of Table 
[T] shows the statistics for the whistler wave amplitudes. The values in each column represent the 
fractional number of events with values satisfying the criteria for each column with at least one 
electric (top row) or magnetic (bottom row) field component. There were 244 whistlers observed 
with at least one electric field component and 65 with at least one magnetic field component. Below 
these ratios are the L-shell ranges over which the whistlers with the corresponding amplitudes were 
observed. Note that a significant fraction, ~40%(~62%), of the electric(magnetic) field whistlers 
have >50 mV/m(>0.5 nT) amplitudes. 

The bottom part of Table [T] shows the wave amplitudes versus AE index statistics as the mean 
± the standard deviation of the mean. The electric(magnetic) field whistlers have large wave am- 
plitudes (>110 mV/m or >1.6 nT) for moderately strong AE (200-400 nT), suggesting substorm 
injections may provide some of the free energy for the observed whistlers. The mean amplitudes 
appear to decrease in the bottom part of Table [TJdue to the limited number of samples for AE > 600 
nT, as seen by the trend in Figure |3] Note that for 10 < AE < 200 nT (see Figure [3]), the average 
wave amplitudes are -0.40 ± 0.04 nT(-31.54 ± 1.39 mV/m), much smaller than for AE > 200 nT. 

Although, much of the literature focuses on the relationship between whistlers and higher en- 
ergy electrons |Li ef a/.||2010| , comparison to electron distributions with energies < 30 keV indicate 
that the properties of the waves vary with the shape of the distributions in this energy range. Also 
previous work has shown that a Cerenkov resonance can excite whistler waves \Kellogg et al. 2010 



Kumagai et al. 1980 Singh 1972 Starodubtsev et al. 1999| , which has much lower energies than 



the typical cyclotron resonance invoked. Note that the effect of Landau damping by electrons in 
this energy range on propagation angles of the whistlers has also been studied \Bortmk et aZ.||200^ . 
Below we will discuss estimates of electron/whistler parallel resonance energies for a range of elec- 
tron pitch angles. The total electron energy is at a minimum for electrons with zero pitch-angle. For 
instance, if a particle has a pitch-angle of 45°, the total particle energy is twice the parallel energy. 
We will discuss those consequences below as well, but first we will discuss just parallel resonance 
energies. Although there are some uncertainties in our density estimates, as explained above, we 
found that most of the whistlers were observed under conditions where they were resonant with 
electrons within the energy range of the EES A detectors. The well known nonrelativistic parallel 
resonance energy [Kennel and Engelmann 1966| of a particle is given by: 



COS0M- — 



CO 

a; 



(2) 



where Bo is the magnitude of the magnetic field, n,. is the ambient plasma density, n„ is the electron 
cyclotron frequency, 0^ is the wave propagation angle with respect to the magnetic field, and m = 
(Landau), —1 (normal cyclotron), or +1 (anomalous cyclotron) for the different resonances. 

For the two examples shown in Figure |2]the necessary parameters needed to calculate resonant 
energies are the following: f - 900-2200 Hz(800-2100 Hz), 0,^ - 73°(39°), N, - 168 cm-3(71 
cm^^), and |B| — 452 nT(168 nT) for the waveform on the left(right). Given the above parameters, 
the Landau resonant energies from the cold plasma dispersion relation are ~ 560-720 eV( 170-240 
eV) and the normal cyclotron resonant energies are ~ 16.2-96.2 keV(0. 37-4.02 keV) for the wave- 
form on the left(right). Of course if our estimates of either the density or propagation angle are 
wrong we receive slightly different results. We know that the PESA High detector saturates in the 
radiation belts and under estimates the densities. Thus, increasing the density lowers both the Lan- 
dau and normal cyclotron resonance energies for both the whistler at 1998-11-13/18:20:59. 590 UT 
and the one at 2000-04-10/03: 10:43.077 UT. Lowering the propagation angle has the same effect, in 
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the cold plasma limit with all other parameters in Equation |2] being constant, since E||„, is roughly 
proportional to 1/cos Ots- 

Particles with a nonzero pitch-angle have a higher total energy than given by Equation|2] This is 
important because the detectors measure the total energy of a particle in a specific angular bin, not 
just the parallel energy. Thus, if the particles have large pitch-angles and Ey^j < 30 keV, it is possible 
that the detector would not measure the particle. This becomes an issue at high pitch-angles (>45°) 
which encompasses a significant fraction of the total distribution when anisotropic. However, the 
waves we observe are highly oblique (as seen below) and thus have large longitudinal electric fields. 
The consequence is that the Landau resonance energies may be more important than the cyclotron 



Tesonance energies \ Kellogg et al. 2010 Kennel and Petscheck 1966 Kumagai et al. 1980 Staro 



dubtsev et al. 1999). Thus, we argue that the distributions produced by the EES A detectors can be 



relevant to the observed whistler waves. 

Recent simulation work by \Sauer and Sydora\ p010[ found that oblique whistler waves can be 
excited by electron beams if the beam speed, V;,, is at least greater than twice the electron Alfven 
speed, = BJy/JiJiJn^, where is the ambient magnetic field magnitude, jUo the permeability 
of free space, n„ the background plasma number density, and m^ the electron mass. However, the 
beams rarely exceed 6000 km/s and the corresponding V^f is typically > 50,000 km/s. Thus, it is 
unlikely that the observed beam/flattop-related whistlers are driven unstable by the mechanism pro- 
posed hy \Sauer and Sydor a |2010|. We also note that the electron distributions used in this study 



have energies almost entirely below the energies expected by ^Sauer and Sydora | 2010) . Thus, a 



beam may exist above 50,000 km/s but we are unable to observe it in the EESA distributions. 



4 Results and Discussion 



We present a statistical study of the properties of large amplitude whistlers and comparison of wave 



amplitudes to AE-index to further characterize the waves described by Kellogg et al. |2010|. We 
have also presented the largest amplitude (>8 nT peak-to-peak) whistler wave measured by a search 
coil in the radiation belts. The highest occurrence probability was within ±6 hours of midnight in 
the radiation belts between 3 < L < 6. In addition, we describe the dependence of whistler charac- 
teristics on the shape of the electron distributions for energies below ^30 keV. However, one should 
note that because the Wind mission was not focused on the magnetosphere, the coverage is limited. 

We also show that very large amplitude whistlers are common in the radiation belts, consistent 
with recent results \Cattell et al. 2008[ Cully et al. 2008 Kellogg et al. 2010 1 using high time res- 



olution field detectors. Many previous studies that focused on time-averaged spectral intensity plots 
underestimated the amplitude of these bursty whistlers while the satellites with waveform capture 
capacities have only recently made measurements in the radiation belts. In addition to their com- 
mon occurrence, simulations have shown that these waves could be capable of producing significant 



2008 Cattell et al. 



changes in the relativistic electron fluxes of the radiation belts | Bortnik et al. 
[20081 . 

The majority of the waveforms observed in our study have f/f„ < 0.5, occurred within L = 6, 
simultaneously with anisotropic electron distribution functions below 30 keV, and the majority of 
the largest amplitude whistlers observed occur during magnetically active periods (AE > 200 nT), 



consistent with recent observations \ Li et al. 2010J . The correlation with AE suggests the source of 
free energy may be due to plasma injections from the geomagnetic tail. Examination of LANL low 
energy (30-300 keV) electron summary data (not shown) showed large injections on the nightside 
at geosynchronous orbit shortly before Wind observed the two whistlers in Figure |2j supporting the 
idea that plasma injections may provide some of the free energy for large amplitude whistler gener- 
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ation. Because the SST detectors were designed for the solar wind, we could not identify features 
in the high energy (>70 keV) electrons associated with the observed whistlers. Regardless of what 
the source of free energy may be for the observed waves, we observe differences in wave properties 
depending on the type of electron distribution functions they are observed with. 

The whistlers associated with anisotropic electron distribution functions in the radiation belts 
exhibited broader frequency peaks than the whistlers associated with beam/flattop distributions in 
the outer magnetosphere. For the 45 TDSS samples in which we were able to determine Oks, we 
found that the waves are often very oblique with Ote's ranging from 0° < Oi^b < 90°, occasion- 
ally more oblique than the whistlers reported by Cattell et al. |2008| and Kellogg et al. | 2010) . 
The anisotropy-related whistlers in both regions were more oblique than the beam/flattop-related 
whistlers observed in the outer magnetosphere. The differences in propagation angle and width of 
frequency peak may be due to differences in the sources of free energy. Some theories suggest that 
whistlers are generated with parallel wave vectors near the magnetic equator and they become more 
oblique as they propagate to higher latitudes. So we examined the relationship between Q^g and 
^asM- We found no correlation between and XasM^ so the waves do not appear to be produced 
as parallel waves which propagate to the satellite and become more oblique as they propagate due 
to dispersive effects. The anisotropy-related whistlers in both the radiation belts and outer magne- 
tosphere are highly oblique, which results in a large longitudinal electric field. When this electric 
field becomes sufficiently large, the typical cyclotron resonances can become less important than a 
Cerenkov resonance \Kellogg et al. 2010[ Kennel and Petscheck 1966 1. Recent studies have found 
that anisotropic bi-Maxwellian electron distributions can produce whistlers with peak growth rates 
for di,B ^ 50°, suggesting that oblique whistlers may be driven by the anisotropic electron distribu- 
tions observed \Hashimoto and Kimura\ [T98 1 [ ^ellogg et al.\ |20 1 0[ [5c/?n ver et al.\ |20 1 0[ . 

Though we observe distinct differences in the whistler characteristics for different electron dis- 
tribution functions, we cannot definitively state that the observed features are the source of free 
energy for these waves. \ Kellogg et a/.|pOTO| presented a warm plasma instability analysis for one 
anisotropic electron distribution function measured by EESA Low near the whistler observed in 
Figure 4 on the right finding the distribution to be unstable. The analysis showed that maximum 
growth occurred at real frequencies corresponding to the observed wave frequency, but the growth 
rate (0.08 s) and propagation angle (^10°) were too low. The instability was found to be a Cerenkov 
resonance \Kennel and Petscheck\ \i966 \ near 150 eV, not the typical cyclotron resonance, due to the 
highly oblique nature of the waves. Owing to the long duration necessary for the EESA Low instru- 
ment to collect a full 4n steradian distribution, one can assume that the instantaneous distributions 
are less isotropic. Thus, \Kellogg et aL\ |2010 1 performed a second instability analysis with a more 
anisotropic hot electron component finding the growth rate (1.5 ms) and propagation angle (up to 
^60°) to be approximately the same as the measured values (they measured ~ 61°). Although, 
the real frequency produced by the analysis did not match the observed frequency. They definitively 
showed that the distribution was unstable to whistler waves through a Cerenkov resonance. Mean- 
ing, the parallel Landau resonant energies estimated by Equation[2]fall well within the energy range 
of the EESA detectors. 



5 Conclusions 

To summarize, we observed 244(65) whistlers inside L = 15 using electric(magnetic) field data from 
the Wind spacecraft; ^40%(^62%) of the waves have peak-to-peak amplitudes of >50 mV/m(>0.5 
nT). The majority of the whistlers had f/f„ < 0.5 and were observed within ±6 hours of midnight 
between 3 < L < 6. The waves had a broad range of propagation angles, 0° < 0^ < 90°, and 0(.g 
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was not correlated with Xgsm- Thus, assuming an equatorial source, the wave propagation angles are 
not a consequence of propagation effects. 

One whistler observed by the Wind search coil had a wave amplitude over two orders of mag- 
nitude above previous observations, large enough to saturate the detector. Also, we only measured 
one component of the magnetic field magnitude for this wave. Thus, the wave amplitude of ^8 nT 
peak-to-peak and our Poynting flux estimate of ^300 /iW/m^ are underestimates of the true values. 
A quasi-linear calculation of the time scales necessary for a wave with this amplitude to accelerate 
plasma sheet electrons to ^1 MeV, assuming 1% efficiency, yields a time scale of roughly 0.5 sec- 
onds, much shorter than previous estimates based on small amplitude whistlers which were on the 
order of days. Using the estimates of maximum change in kinetic energy fTom \Omura et al. (2007), 
we found the 8 nT whistler to be capable of producing ^61 MeV electrons through relativistic turn- 
ing acceleration. 

Although we were not able to determine a source of free energy for the observed waves there 
were physical differences in the whistlers associated with different electron distributions. In addi- 
tion to the anisotropy-related whistlers having larger Af/f than the rest of the whistlers, we observed 
differences in wave characteristics with the following electron parameters: T±,eJT'\iel^ Timh^Tieh^ 
Tel, and T^^ (where EL = EESA Low and EH = EESA High). The whistlers observed in the radi- 
ation belts simultaneously with anisotropic electron distributions had larger Af/f than the rest of the 
whistlers observed. We also examined the LANL low energy (30-300 keV) electron summary data 
and found that our largest events were associated with large substorm injections, consistent with the 
increase in wave amplitudes with increasing AE-index. Our study adds to the mounting evidence 
that very large amplitude whistler waves are an important phenomena in radiation belt dynamics. 
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Table 1: Wind Petal Orbit Whistler Statistics 



Whistler Location and Particle Statistics (L < 15) 


Electron Distribution Function 


Inside of 


Between Magnetopause 




Type 




Radiation Belts 


and Radiation Belts 




Anisotropic 




172 


42 




Beam/Flattop 




2 


28 


Whistler Amplitude Statistics (L < 15) 


Whistler Waves with |Eh | 


>100mV/m 


>50 mV/m 


>30 mV/m 


#/(Total # measured with |E,v |) 


48/244 


97/244 


129/244 


Range of L-Shells 


3.8-5.8 


3.7-6.3 


3.7-7.9 


Whistler Waves with |B„ | 


>1.0 nT 


>0.5 nT 


>0.3 nT 


#/(Total # measured with B,^ |) 


17/65 


40/65 


44/65 


Range of L- Shells 


4.0-8.8 


4.0-11.3 


4.0-11.3 


Whistler AE (nT) Index Statistics (L < 15) 


Type 


200 < AE < 400 


400 < AE < 600 


600 < AE < 800 


800 < AE < 1000 


|E, 1 (mV/m) 


110.88 ±6.45 


100.56 ± 5.69 


96.43 ± 4.75 


40.35 ± 1.86 


|B. 1 (nT) 


1.61 ± 0.27 


1.09 ± 0.04 


1.65 ± 0.03 


0.92 ± 0.03 


#/(Total#of |E„, 1) 


34/217 


60/217 


14/217 


35/217 


#/(Total#of B„ ) 


13/52 


9/52 


2/52 


13/52 
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11-13-1 998/1 3:00:00 -11-14-1 998/02:00:00 UT 



5000 




10 b 

Uni.Time 14:00 
MLT 17.3 
L-Shell 9.80 
Rad. Dist. 8.10 



27 keV 
40 keV 

bt) KeV 
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520 keV 



71 keV 
128keV 
196 keV 
ii3 keV 
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1017 keV 
2073 keV 
4439 keV 
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22:00 
7.10 
9.40 
9.40 



00:00 
7.90 
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02:00 
8.40 
14.9 
14.9 



Figure 1: Example of a Wind petal orbit for the time range of 1998-11-13/13:00:00 UT through 
1998-11-14/02:00:00 UT. The top panel is a time averaged spectral intensity plot from the 
WAVES/TNR receiver with the local f„ (black line) and f„/, (white line) over plotted. The second 
panel is the GSM components of the magnetic field and magnitude. The third and fourth panels are 
omni-directional number fluxes of high energy electrons and protons, respectively, from the Wind 
SST instruments. At the bottom of the plot are tick mark labels of MLT, L-Shell, and radial distance 
along with the local UT. 
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TDSS: 11-13-1 998/1 8:20:59.590 UT TDSS: 04-1 0-2000/03:1 0:43.077 UT 




40 80 120 160 200 240 40 80 120 160 200 240 



Time (ms) Time (ms) 

Figure 2: Examples of two whistler waves, in instrument coordinates, observed on the perigee passes 
of Wind on 1998-11-13 and 2000-04-10. The examples were taken from the TDSS instrument 
and filtered (frequency range shown in top panels) to remove superposed low and high frequency 
signals. The top row shows the Y-component of the magnetic field (blue), second row shows the X- 
component of the electric field (red), third row shows the Y-component of the electric field (green), 
and the fourth row shows the Z-component of the electric field (magenta). The amplitudes of each 
component are marked by the vertical black arrows. The TDSS samples for these two events were 
taken at 7.5 kHz. 
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Peak Whistler Amplitude vs. AE-lndex 




200 



400 600 800 
AE-lndex [1 min, nT] 



1000 1200 



Figure 3: Peak whistler wave amplitudes versus one minute AE-Indices for the whistlers. The top 
panel corresponds to the peak whistler amplitudes measured with the Wind search coil (ranging from 
~0.01-2.0 nT) and the bottom panel is the peak whistler amplitudes measured with the Wind electric 
field antennas (ranging from ^0. 1-350 mV/m). The range of both horizontal axes is 10-1200 nT for 
the AE-Indices. 
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Velocity (1000 km/s) 

Figure 4: Example cuts of electron distribution functions for beam/flattop (Type A), beam (Type B), 
and anisotropic (Type C). The parallel(perpendicular) cuts of the electron distribution functions are 
seen as the red(blue) solid(dashed) lines. The green line corresponds to the one-count level for that 
particular distribution. 
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Y-GSM (Rg) vs. X-GSM (Rg) 
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Figure 5: The top four plots show 12 Wind petal orbits over layed and color coded (labeled in each 
panel) with whistler wave locations (defined by color coded *) projected onto the XY-GSM plane. 
For each petal orbit, the start(end) point is identified by a A(<^). Below the four petal orbit plots are 
two histograms of GSM-Latitude, )iasM (degrees), for the 239 TDSF and TDSS samples on the left 
while on the right is the histogram of (degrees) for the 45 TDSS used in MVA. Each histogram 
panel contains three histograms separating the whistler waves by the concurrent electron distribu- 
tion function observations and magnetospheric location. The plots are organized in the following 
order: (top panel) beam/flattop electron distribution functions in the outer magnetosphere, (middle 
panel) anisotropic electron distribution functions in the outer magnetosphere, and (bottom panel) 
anisotropic electron distribution functions in the radiation belts. The vertical axis is the percentage 
of whistlers in each event category. 
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